Spectroscopy and 3D imaging of the Crab nebula 
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ABSTRACT 

Spectroscopy of the Crab nebula along different slit directions reveals the 3 dimensional 
structure of the optical nebula. On the basis of the linear radial expansion result first 
discovered by Trimble (1968), we make a 3D model of the optical emission. Results from 
a limited number of slit directions suggest that optical lines originate from a complicated 
array of wisps that are located in a rather thin shell, pierced by a jet. The jet is certainly 
not prominent in optical emission lines, but the direction of the piercing is consistent with 
the direction of the X-ray and radio jet. The shell's effective radius is ~ 79 seconds of 
arc, its thickness about a third of the radius and it is moving out with an average velocity 
1160 kms"^. 

Subject headings: ISM: individual (Crab Nebula) — ISM: kinematics and dynamics — supernova 
remnants — techniques: spectroscopic 



1. INTRODUCTION 

The Crab nebula was an object of intense interest even before its pulsar was discovered in radio by 
Staelin and Reifenstein (1968) and in the optical by Cocke, Disney, and Taylor (1969). Trimble (1968) 
first measured proper motions of nebular components and, combining them with spectroscopic data 
taken by Guido Miinch (1958), discovered that they all point to a very small, almost point-like origin 
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in the past. Velusamy, Roshi, & Venugopal (1992) obtained a similar result from radio images of the 
nebula. Both observations are consistent with an outward pointing velocity field with the proper motion 
magnitude proportional to the distance from the origin. Trimble (1968), Wyckoff and Murray (1977), 
Caraveo and Mignani (1999) pointed out that the position of this origin is displaced with respect to 
that of the pulsar. Furthermore the proper motions of nebular components are about 10% faster than 
expected if the nebular debris were expanding freely. In radio images the linear expansion is not as 
obvious and shows some anomalies further away from the pulsar. More recently, attention has focused 
to phenomena connected to Crab's high-energy emission, its jet and to the very rapid expansion of 
inner wisps discovered optically by Hester et al. (1995) and in radio by Bietenholz et al. (2001). It 
has become increasingly clear that the Crab pulsar and its nebula are connected through a multi layer 
structure of different predominant energy ranges. To understand the mutual interrelationship of these 
layers, it is useful to have a three dimensional picture of at least part of the structure. 

The first attempts to understand the 3D structure of the Crab nebula go back to Chevalier and Gull 
(1975) who studied the nebula in the light of all 13 important visible emission lines by using narrow band 
filter imaging. They found a structure "consisting of dense filaments from which thin clumpy sheets of 
gas fan out radially". Clark et al. (1983) studied spectroscopically the radial velocity field across the 
nebula to investigate the nebula 3D properties and find the line emitting region to be a thick hollow 
shell composed of bright inner and faint outer components. They also see the filaments as generally 
circumferential, but radial "spokes", connecting the inner and outer shells. They conclude that the 
nebular synchrotron emission is confined within the inner brighter shell. A comprehensive attempt to 
understand the geometry, composition and mass of the Crab nebula has been done by MacAlpine et 
al. (1989) and (1996). They combine interference filter imaging and long-slit spectroscopy to obtain 
information on composition and distribution of gas in the nebula. On the basis of a few selected spectra 
and using model calculations MacAlpine et al. (1989) deduce that coUisional contributions to observed 
exceptionally strong helium lines are insignificant. Using further model assumptions they estimate the 
combined mass in the filaments to be in the range 6-9 Mq. On the basis of interference filter and Fabry- 
Perot imaging Lawrence et al. (1995) interpret the [OIII]A5007 emission to have a spatial relationship 
to a " high-helium torus" . 

The possibility to make a 3D representation of the line emitting regions in the Crab nebula is offered 
by the proper motion result of Trimble and others. On the basis of this result one is led to speculate 
that all the three components of the velocity vector field describing the motion of nebular components 
point in the outward radial direction since the motion is ballistic. According to this assumption the 
speed of this motion is proportional to the distance from the origin, i.e. we assume: 

v = r/T , (1) 

where v is the velocity vector at the position r measured from the center of the nebula, and T is time 
elapsed. If the ballistic assumption was fully satisfied, T would be the time since the explosion in 1054, 
but, as Trimble's proper motion results show the possibility of some acceleration of the debris in the 
past, we consider T as a model parameter that may differ somewhat from the time since the explosion. 
Observing spectroscopically a small patch of the nebula at the angular position {$^, rf) from the pulsar, 
one expects to find multiple Doppler shifted emission lines, each Doppler shift corresponding to a 
different radial velocity component (f^) of the cloud. According to (1), such a component corresponds 
to a line of sight coordinate z — {vz — Vzc)T, while the other two components of the vector r are 
X = {$,—$,c)D and y = {rj — ric)D, where Vzc and (^^c, ■Cc) are the velocity and (angular) coordinates of the 
center of the cloud, while D is the distance to the nebula. Locating line emitting regions in this way. 
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one can construct a tentative 3D image of (line emitting) nebular components. The only parameter 
critical for the shape of such a reconstructed image is the ratio T/D. We choose to consider it as a free 
parameter and construct 3D models of the nebula for different values of it. We find it remarkable that 
it is possible to find a value for this parameter such that the 3D model of the nebula in the radiation 
from the strongest optical emission lines ([Nil], [SII], and H^) appears as embedded inside a relatively 
thin, pierced spherical shell with a well defined inner radius. We consider this particular ratio T/D 
of some interest, yet, we would like to emphasize that its implications for the distance to the nebula 
depend strongly on the ballistic assumption. 



2. DATA 

The spectra of the Crab nebula were obtained on the 18*'* of January 2001 (slit orientation 0° and 
90°, high resolution), on the 4*^^ (0° and 90°, low resolution ), S*'* (18°, 36°, 54° and 72°, low resolution), 

(108°, 126°, 144° and 162°, low resolution), (36°, 72°, 108° and 144°, high resolution) and 8*'' 
(18°, 54°, 126° and 162°, high resolution) of January 2002 with the 2.12 m telescope of the Observatorio 
Astroffsico Guillermo Haro in Cananea, Mexico, equipped with the Boiler and Chivens spectrograph. 
All observations were done under photometric conditions. Spectra were taken in a low resolution mode 
of 1.6A per pixel and a 250/xm (a; 4 pixels ~ 2 arc seconds) wide slit, and in a high resolution mode 
of O.43A per pixel and a 50//m slit. The spatial FWHM of the pulsar profile was varying between 2 
and 3 arc seconds during observations and was thus wider than the projected width of the sht. The 
slit length was 560 pixelsxO.46 arc seconds per pixel in both cases, however, all the pixels were not 
exposed properly (see Fig. 1). In the first case the effective Doppler width of the slit was 300 km s^^ 
and in the second it was about 16kms~^. Five spectra (20 min exposure each) were taken with the low 
resolution and three (30 min exposure each) with the high resolution for each of the 10 slit orientations 
(■0s = 0°, 18°, 36° . . . 162° with respect to the EW orientation). The low resolution spectra cover 
the spectral region between 5400A and 7000A and give an average signal to noise ratio (S/N) of more 
than 200 per pixel. The high resolution mode covers the region between 6400A and 6840A and yields 
an average S/N of only about 2. Thus, in the low resolution mode the nebular continuum is clearly 
detected with the S/N of up to about 40 per pixel, while in the high resolution it is just below the 
detection threshold. As an example, reduced (bias subtraction, sky ffat, median average) high and low 
resolution spectra for the slit orientation 36° are presented in Fig. 1. In Fig. 2 we also show a rescaled 
superposition of the two, where the high resolution spectrum is colored in red and the low resolution 
one in green. A trace of the spectrum (at about 75 arc sec from the pulsar, below in Fig. lb) is shown 
in Fig. 3. The rest frame wavelengths of a number of known nebular spectral lines arc marked in gray 
and the Doppler shifted lines of nitrogen, sulfur, hydrogen, oxygen, and helium are marked in blue 
(— 1047kms~^) and red (+1729 kms^^). The scale on the ordinate axis is in photon counts per pixel 
per 20 minutes. Finally, Fig. 4. shows the continuum trace along the V's = 36° slit at approximately 
6OOOA, i.e. in the fiat region between the Hel and [01] lines. The dot shows the position of the pulsar, 
where the continuum value interpolated along the dotted line is calibrated with respect to previous 
measurement (Carramiiiana et al. 2000). Some features in this and other continuum traces can be 
recognized as belonging to stars that can be identified in the image of the nebula (the peak just next 
to the right of the pulsar belongs to the star k, 6 arc seconds away). We use this information to secure 
accurate positions of the slits with respect to the nebula. The positions thus derived are shown in Fig. 5. 
The (small) corrections with respect to desired positions, obtained by these data have been taken into 
account in further data analysis of low resolution spectra. This procedure could not be apphed for the 
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high resolution spectra because of the low S/N. In this case we assume that the slits crossed the pulsar 
accurately and that the transparency of the sky was constant. We expect that the respective error 
in position and in flux is comparable to the shifts and calibration factors obtained in low resolution 
spectra (±4 arc seconds in position and about ±10% in flux). In constructing images we use the more 
precise position data from low resolution spectra and use high resolution spectra to reduce the velocity 
calibration errors of the low resolution spectra. 

An analysis of all the spectra conflrms the impression apparent from Fig. 3, that the intense [Nil] 
and [SII] doublets, Hq, line and also the weak [01] doublet always occur at a common red-shift, relative 
intensities varying across the nebula. The forbidden [NII]5755 line is almost always clearly absent, but 
in some instances its presence may be masked by the double atmospheric feature just to the left and to 
the right of the red [Nil] 5755 mark (compare also Fig. 1). Only a close scrutiny revealed the presence 
of this weak hne in the 36° spectrum at the red-shift of the strongest [Nil] line. It is quite difficult to 
automatically obtain good estimates for the ratio (/[6548] -|- /[6583])//[5755] from these spectra, since 
the continuum is not well defined in this region because of the mentioned presence of weak atmospheric 
lines. Checking spectral traces by hand and judging the intensity of the /[5755] line by the value at the 
pixel corresponding to the velocity of the /[6548] and /[6583] line, we estimate some values for this ratio 
between 185±35 (7[6548] +/[6583] = 259photon counts/pixel), 210±35 (7[6548] +7[6583] = 487photon 
counts/pixel) and 310 ± 100 (/[6548] + /[6583] = 193photon counts/pixel) for traces, where the J[5755] 
line is visible. This value is consistent, even somewhat lower than that found by Davidson et al. (1982). 

Among the nebular lines observed in our spectra, we take only the five most intense as characterizing 
optical emitting regions: [Nil] {i = 1, 2), Hq, {i = 3) and [SII] {i = 4, 5). To determine 
their common Doppler shift and strength we proceed as follows: for each slit orientation we scan the 
spectra along all slit length positions m (10 pixels < m < 537 pixels for high resolution spectra and 
43 pixels < m < 545 pixels for low resolution spectra; the slight misalignment of the slit with respect 
to the CCD axis is taken into account) and obtain raw spectral traces S'r(A, m|-?/'s). These arc fiux 
calibrated by making use of the fact that all slit positions cross the pulsar. Thus we set the 6000 A 
interpolated nebular continuum at the position of the pulsar to the same reference value for all spectra 
and analyze all spectra in units of this reference value (see Fig. 4)^. Pure hne spectra Sl{X,m\^jJs) are 
then obtained by subtracting the continuum. It is automatically defined for each trace (m) as the 

best n*^ (n=10) degree polynomial fitting the points {Xf'\ Sr{Xf'\m\'ilJs)), where A-'^'* are determined by 
inspection of the whole spectrum at the wavelengths free of spectral lines^. Lines and their common 
Doppler shifts are finally identified defining a Doppler velocity dependent weight for each of the five 
lines in the spectrum as: 



and find on the interval —2000 km s ^ < f < 2000 km s ^ the velocity Vmax at which Q{v) is maximal. 
Such a velocity weight function is large only if both lines of doublets are present and it has been found 



^We estimate that the relative error of this recahbration is of the order of ±3%, about the size of the dot in Fig. 4. 

^The recalibration of the spectrum and continuum subtraction is possible only for the low resolution spectra where the 
continuum is observed with sufficient signal to noise ratio. In high resolution no continuum is observed and no subtraction 
is applied, therefore, no common calibration intensity is present. 
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Choosing the appropriate value for the parameter a we form the velocity weight function: 

n(f) = Wi{v)w2{v) + WcW3{v) + W4{v)w5{v) 



(3) 
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experimentally on our data sample (some 55,000 spectral traces) that its maximum gives a velocity 
that can be considered the common velocity of the three components (H, N and S). The characteristic 
strength Wc introduced with the Hq, line is chosen so that the contribution of this line to the velocity 
weight function is in the average comparable to the contribution of the two doublets. The Gaussian 
lines, Doppler shifted to Vmax arc then subtracted from Si{X,m) and the procedure is repeated on 
the remaining spectrum until Q{v) no longer shows significant maxima. The important parameter a, 
chosen as constant in the reduction of a given spectrum, is selected before the automatic procedure 
starts. This is done through trial and error until most randomly chosen line multiplets subtract best 
in a single subtraction. In the low resolution case the best a value is found to be essentially the 
effective sht width. This is expected, since no hues show evidence of being resolved. In this case the 
subtractions of all the five lines belonging to a common velocity multiplet usually leave only little more 
than the expected photon counting statistical noise, except in few cases where different Doppler shifted 
components of Hq, and [Nil] overlap. The weights Wi{vmax-,^) form what we call reduced spectra and 
the set of the five weights for a given Vmax and m is called here a velocity component. Three examples 
of successive subtractions are shown in Fig. 6. In this way we detected 43843 velocity components in 
all the 10 low resolution spectra, in the average 8.5 components per spectral trace. A great majority of 
these are of no consequence, since their intensity is negligible. 80% of light is contributed by only 5300 
(12%) of different velocity components, and 20% account for 90% of light. In other words, 8770 velocity 
components (1.7 per spectral trace) account for 90% of hght. The remaining 80% of detecte velocity 
components contribute little, so they are barely seen in the reconstruction. Some of them belong to line 
wings that are not subtracted during the first run, to weak independent lines, to residues of imperfect 
subtraction which are due to either photon counting noise or small errors in wavelength calibration, but 
also due to the fact that different lines sometimes overlap and their components become confused. It is 
difficult to precisely judge the relative contribution of these effects. We beheve they are best illustrated 
in Fig. 6, which shows typical examples for reasons of imperfect subtraction to the level of a few percent 
and in Fig. 7, where the observed spectrum and its reconstruction from detected velocity components 
can be compared. 

The high resolution spectra clearly show lines of different resolved widths. A drastic example of 
both an unresolved [SII] doublet at — 115kms~^ and a 150kms~^ broad [SII] doublet at 814kms~^ is 
found in the spectral trace along the hne in Fig. la and is shown in Fig. 8^. In this case it would be 
advantageous to take cr as a free parameter for each system of lines. However, since most lines are 
neither as narrow nor as wide as the example shown here and the S /N is quite low, we usually choose a 
moderate a = 40kms~^ which subtracts most lines rather well. Only the broadest and strongest lines 
produce additional sidebands next to the main peak in reconstruction. A clear example of such a broad 
hne is clearly seen in Fig. 10 on top of the t/j — 162° and also on top of the ■0 = 0° spectrum. 

Spectra reconstructed from velocity components are shown in Figs. 9 and 10, where the intensity 
of the weights is color coded as indicated by color triangles in the middle. The intensity is defined as 
I— J — ^-TT\{wi + + ^3 + W4 + uis), where Imaxii^) is the maximum intensity for the given reduced 
spectrum. In units of interpolated continuum intensity at the pulsar position Imaxi'^P) are: 42.5, 22.4, 
61.7, 41.7, 38.3, 46.0, 30.1, 43.0, 57.5, 74.0 in order from = 0° to ^ = 162°. The colors are assigned 
with respect to the "average bright spot mixture of intensities", which is coded gray (center of the 



^Note that the Doppler width of thermal broadening for hydrogen at the temperature 8000K is 11.5 kms~^, so our 
spectral resolution is not sufficient to show line broadening due to this effect. We expect that at this level of resolution 
the width of the line tells mainly about the spatial extent of the radiating gas. 
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triangle, {H : N : S}={0.18 : 0.42 : 0.40}), while the edges of the triangles correspond to pure H, 
N or S as indicated on top in Fig. 9. Comparing Figs. 9 and 10 (compare also Fig. 1) one notes 
the advantages and disadvantages of high and low resolution. The dominating features arc clear and 
accurately correspond in both figures. However, the global features of the radiating gas distribution are 
better appreciated in low resolution due to higher signal to noise ratio. On the other hand, only the 
high resolution shows that the radiating gas is organized in exquisitely fine sheets, tubes or blobs. 

Reduced spectral intensities are considered the line intensities at the position: 

^ = mcos'^ + A^('0), 77 = msin'^ + A77('0), v — v^ax ^'^{i^) (4) 

of the 3D image and are stored accordingly. Here m — [m — mo{il})] x 0.46 arc sec/pixel and mo{ip) is 
the value of m in the spectrum 5'r(A, m) where the spectral trace belonging to pulsar light is strongest. 

and Ar]{ilj) are coordinates of the vector connecting the slit position (mo('?/')) with the center 
of the pulsar in the image; see the lower part of Fig. 5 (for A,^(0) and Ari{0)). The A^, Arj values 
are of the order of 2 arc seconds, except for the slits ip = 72° and 162°, where the star next to the 
pulsar was mistaken for the pulsar. The velocity calibration error of low resolution spectra Av{ip) was 
determined by graphically correlating low resolution and high resolution synthetic spectra, assuming 
that high resolution spectra velocity error is negligible. The values found are between —180 km s^^ and 
Okms~^. From high resolution spectra A,^(^) and Af]{ijj) and Av{ip) can not be obtained, but the 
comparison with low resolution spectra suggests that the pointing error was not larger than the one 
reached in taking low resolution data. 



3. IMAGES 

On the basis of the spectral information alone we first reconstruct a) the continuum image of 
the nebula at 6000 A and b) the image of the nebula in the light of the five strong lines. They 
are shown together with the optical, radio and X-ray image in Fig. 12. The black and white image 
corresponding to the 6000 A continuum (bottom left) is formed by assigning the values for the continuum 
intensity at measured points {^k, Vk) and values in between are assigned by linear interpolation, i.e. if 
^^(A, Am + mQ{i/js)\'^s) = Si and Sr{\, Am + mo{i/js)\'^s + Ail!) — S2, then Sr{\, Am + mo{ilJs)\ips + S'lp) = 
Si + {S2 — Si)-^, where A^ = 18° and < Sip < Aip, and the coordinates ^ and rj belonging to points 
m — Am + mo{ip), ip = ipg + Sip are calculated according to eq. 4. 

The emission lines image (bottom right) is formed in a similar fashion, but in color so that the 
green intensity is proportional to the intensity of nitrogen regardless of velocity, the red component 
comes from Hq, and the blue from sulfur. The color contributions are normalized so that the strongest 
line of any color gives the maximum display value for this color (255). 

In order to appreciate the level of detail that such synthetic images can give, we also display in 
Fig. 12 a blend of the visual image and the X-ray image (top left), the radio image (top center) and 
an enhanced visual image (top right) of the nebula. We find the continuum image reminiscent of the 
X-ray image, while the N, H, S image clearly highlights the filament structure embedded in the nebula, 

as does the radio image. 

The 3D image can be understood on a computer screen if the structure is rotated and displayed as 
a movie (which can be obtained in the electronic version of this paper; see also Fig. 13) which is based 
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on low resolution data. The color code is the inverse^ of that used in Fig. 9 and the intensity is coded as 
density of radiating fog multiplied by sin(^), where 6 is the angle between the Z axis and the direction 
from the coordinate system origin to the observed cloud. This is to compensate for the difference in 
sampling density. The axes X and Y are oriented as in Fig. 5 and their length is 90 arc sec, while the 
Z axis points away from the observer and its length is 1300 km/s. When observing this animation, it 
becomes quite apparent that the structure is shell hke and that it is possible to adjust the T/ D value 
so that it looks round when rolled on the screen. Therefore, we fit a spherical shell to the distribution 
of line radiators by minimizing the action: 

A[i^, r)„ Vc, A, a] = ^J— J] U{ik - Cc)' + {Vk - Vc)^ + a^{vk - Vc)^ - aV Uk . (5) 

Here Uk is the weight of the k-th radiating blob formed from line intensities^ Wi{^k,Vk,Vk) and the 
variational parameters {Cc,Vc), Vc, A, a are the position of the center of the shell with respect to 
the pulsar, the velocity of the shell center with respect to the observer, the effective shell radius 
and a — T/D respectively. In the case of the low resolution the best fit was obtained for the fol- 
lowing parameter values: {^c: Vc: A, Vc) — (10.4 arc sec, 17.8 arc sec, 79.1 arc sec, 198kms~^) with the ac- 
tion value of (13.2 arcsec)^. In the high resolution case the parameter values are: {^c,Vc, A.,Vc) = 
(11.3 arc sec, 13.8 arc sec, 74.5 arc sec, 189kms~^) for the action vahic of (10.1 arcsec)^. The best fit pa- 
rameters were used to plot projections of the shell on the v — p planes in Figs. 9 and 10. The dynamical 
center of the nebula as found by our fit {^c, Vc, ^c) is plotted as a cross in Fig. 12; both the fine 
resolution and the low resolution spectra obtain center within the box defined by the size of the cross 
(about 10 arc seconds). The best value for a is 0.068 arc sec/ (km s^^). This position of the center of the 
nebula is puzzling. The vector from our center of the nebula to the pulsar points almost exactly opposite 
to the proper motion vector multiplied by time since explosion found by Caraveo and Mignani (1999). 
This could be explained if the nebula was moving with respect to us with twice the velocity of the 
pulsar, which seems unlikely since Crab is considered a galactic disk object. Therefore, we checked the 
assumption that the center of the nebula is where Caraveo Mignani's proper motion result would put 
it, assuming that the proper motion of the nebula is negligible. We constructed a movie rotating about 
this center. In this case wobbling of the 3D image is beyond any doubt. Note, however, that due to slit 
length limitations the radial distribution is cut at ~ 120 arc seconds from the center, where the intensity 
is not completely negligible. This should be taken as a warning that the obtained fit values may not yet 
have a very simple physical meaning. If we assume in (1) that T is the time elapsed since the supernova 
explosion, we obtain 3.1 kpc as the best value for the distance, somewhat larger than the accepted value. 
Finally, in Fig. 14 we show the radial distribution of intensity (proportional to the power radiated in 
a line in the shell between p and p + Ap) in the three line sets as a function of the "radial distance" 
p = \/ — ^c)^ + {v ~ Vc)'^ + Q;^('y — Vc)'^. The sharp rise at a; 70 arc seconds is followed by gradual and 
regular fall off. 



^cyan— >red, yellow— >blue, magenta— >green 

^For high resolution spectra = X)f=i ''^i(S,k,ilk,Vk)- In the case of low resolution spectra the inspection of Figs. 9 
and 12 makes obvious that the p — v plane of the nebula has not been scanned exhaustively, some waining filaments having 
clearly escaped detection since the slit of the spectroscope is not long enough to reach to the farthest radiating regions. 
We note, however (Fig. 9), that the intensity generally decreases with the distance from the shell, so that weighting 
brighter clouds more heavily than the dimmer ones might not displace the center of the nebula, but it does overcome the 
geometric selection problem. Therefore, we cither exclude all clouds with less than a few percent of maximum intensity 
or by weighting Uk proportional to w"^. Both methods give similar results which differ only by about 1% in A, by about 
5 arc seconds in {^cVcCtVc) but almost 30% in a and agree with those from the high resolution spectra. 
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Our 3D image suggests that the volume within ^ 55 arc seconds of the center of the nebula is almost 
void of line emitters. A noteworthy exception is a very weak trail of spots on the line connecting the 
upper left and lower right brightest spot and going through the center of the nebula in the spectra of 
slit orientation 144°, which is the orientation of the X-ray and radio jet in the image (sec Fig. 5). These 
spots are less than 5% of maximum intensity and are likely along the jet. A saturated image of this 
reduced spectrum is shown in Fig. 11. 



4. CONCLUSION 

We obtained multiple long slit spectra of the Crab nebula in the light of the optical emission lines 
of Hq, [Nil] A6548, 6583, and [SII]A6716,6730. The simple assumption of baUistic expansion, expressed 
by (1) and suggested by the work of Trimble (1968) allows a 3D representation of the nebula, although 
through an arbitrary scale in the line of sight. Such a 3D representation suggests that line emitting 
regions are organized in a relatively thin shell, regardless of the line of sight scale. Therefore, we make 
the assumption that the shell is likely to be quasi spherical and adjust the unknown scale parameter 
accordingly. The resulting 3D image indeed appears to fill a spherical shell with an inner radius of 
f» 55 arc seconds, as already shown by results of Lawrence et al. (1995). The shell is pierced by the jet 
in the SE direction and at about 75° with respect to the line of sight. The center of this shell is well 
aligned with the two-dimensional projection of the center of the optical and X-ray continuum emissions. 
It is clearly displaced from the pulsar but appears to be connected to it through the north-western 
extension of the jet. The 3D distribution of the optical continuum and the X-ray emission can not 
be determined, but the observed 2D projections of these are very similar and consistent with the idea 
that the optical and X-ray continuum come from the same, rather fiat, oblate structure centered in 
the spherical shell. The optical continuum emission drops regularly with distance from the center of 
the nebula (Fig. 4), becoming weak when reaching the spherical shell. Based on the absence of the 
[NI1]5755 line, we estimate that the effective gas temperature in the shell is below 7500K. The filling 
of the 55 arc seconds shell seems patchy. We note that some bright spots are only arc seconds wide (cf. 
Fig. 1) and possibly less than 16kms~^ "deep" (Fig. 8), so that we expect that a more detailed high 
resolution 3D view would show that the shell is woven with delicate threads. Note that the small size of 
line emitting regions indicates high emission coefficients in the bright regions; in particular, we estimate 
that the region shown in Figs. 1 and 8 is only about 0.02 ly thick which leads to the emission 
coefficient of the order ja^p 10~^^ ergcm"^ sec~^, which is an order of magnitude higher than the 
average emission coefficient at these wavelengths of the whole nebula including the synchrotron cloud. 
Hints of threads are suggested in the superimposition of the radio image by Vclusamy et al. (1992) on 
the optical image and are clearly seen in the newer radio image by Bietenholz et al. (2001) shown in 
Fig. 12. Note that optical filaments closely map the smoother radio loops. The jet, which is so prominent 
in X-rays and in radio, is marked also by optical continuum emission and seems to be threaded also by 
weak line emission (see Fig. 9, in particular at the 144° orientation and Fig. 11). Line emission from 
regions where the jet hits the 55 arc seconds shell is particularly strong - the upper left brightest spot 
in the reduced spectrum ips = 144° and also i/jg = 162° of Fig. 9. 

In concluding, we would like to call attention to two observations: The roundness and average 
regularity of the spherical shell as expressed by the sharpness and simplicity of the radial distribution 
of line emitting regions in Fig. 14 is quite remarkable. The delicate radiating tubes or blobs filling 
this relatively thin shell closely follow radio loops on large scale, even if, as emphasized by Sankrit et 
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al. (1998) and Hester et al. (1996), on short scale they are dominated by the magnetic Rayleigh- Taylor 
instabilities. The 3D image, even if it still lacks in detail, suggests that the filaments are wound on the 
shell along long smooth threads, possibly following the magnetic field fines. 

This work was done as a collaboration of the Instituto Nacional de Astrofisica, Optica y Electronica 
(INAOE) and the University of Ljubljana and was supported in part by CONACyT grant 25539E and 
by the Ministry of Science, Education and Sport of the Republic of Slovenia under grant number 1554- 
501. The essential support to this project by the technical staff of the OAGH in Cananea is appreciated. 
The movie production is a courtesy of ARTREBEL9. 
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Fig. 1. — (a) High and (b) low resolution of the reduced spectrum for the slit angle 36°. Three 30 
minute exposures for the high resolution spectrum were obtained on 2002 January 8 and five 20 minute 
for the low resolution spectra were obtained on 2002 January 6. High resolution spectrum covers the 
region between 6406A and 6842A (0.43A/pixel) and the low resolution spectrum covers wavelengths 
between 5387A and 7054A (1.6A/pixel). The vertical scale is 560 pixels x 0.46 arc seconds per pixel in 
both cases, but only pixels between 10 < m < 537 arc properly exposed in the high resolution spectrum 
and between 43 < m < 545 for the low resolution spectrum. Continuum is calibrated in the brightened 
vertical region centered at 6OOOA shown in the lower half of the figure (b) . 



Fig. 2. — Low (green) and high (red) resolution spectra in the spectral region between 6415 and 6814 A 
for the slit angle 36° in the vicinity of the [NII]A6548,6583 doublet, H«A6563, and [SII]A6716,6730 
doublet suitably rescaled and superimposed. 



Fig. 3. — Trace Sr{X,m) of the spectrum m = 440 pixels (along the line in Fig. lb). Some important 
lines are marked; the rest positions with respect to local lamp calibration are shown by gray vertical 
lines and parallel blue and red lines indicate blue/red shifted positions by —1047 km s~^ and 1729 km s~^ 
respectively. Absolute velocity calibrations in low resolution spectra are not better than ^ lOOkms"^, 
but comparisons like the one in Fig. 2 show that calibration shifts are the same, to within the resolution 
of the high resolution spectra, in the whole region of [Nil] and [SII] hues of interest. The scale on the 
ordinate axis is in photon counts per pixel per 20 minutes. 



Fig. 4. — Continuum intensity at 6OOOA along the 36° sht (see Fig. lb). The abscissa ticks 

are at 10 arc seconds spacing, and the ordinate shows the intensity as calibrated according to 
(Carramifiana ct al. 2000; Vcron-Cctty ct al. 1993). 1 x 10^^'"^crg scc^-*^ cm^^ A^^ arcscc^^ correspond 
to approximately 10 photon counts per pixel per 20 minutes. The dotted line indicates the continuum 
interpolation and the dot is placed at the intensity corresponding to the interpolated nebular continuum 
at the position of the pulsar, where calibration is applied. 



Fig. 5. — Slit positions with respect to the nebula and their reconstruction. The slit orientation is as 
read from the dial on the spectroscope. The value of m at the pulsar trace in the spectrum is taken 
as reference. The coordinate axis X (^) points to the West (right) and Y (?]) points to the North (up). 
Star traces in the spectrum are noted and the corresponding distances to the reference trace are marked 
as points along a line with a given orientation. Such points (as for example the three black dots in the 
bottom image for slit orientation 0°) are parallelly transported until they correspond with some stellar 
feature in the image. The background image is a combination of three V, R, I filter images obtained at 
the Guillermo Haro observatory in Cananea by Cadez et al. (2001) in the stroboscopic mode to enhance 
the brightness of the pulsar. 
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Fig. 6. — Two successive subtractions of Gaussian multiplets at remaining maximums of the velocity 
weight function. In the three parts of the figure exemplary sections of low resolution spectra around 
the B.a, [Nil] and [SII] lines are shown as border lines between gray and white. The first subtraction 
leaves the red spectral trace, and the black line marks what remains after the second subtraction. 
The two right hand side [Nil], Hq,, and [SII] quintuplets are clearly separated. Noise in the spectrum 
remaining after the first subtraction is roughly consistent with photon counting noise. The remaining 
weak multiplet (red) is again subtracted out to a line whose noise is consistent with the photon counting 
noise. The middle spectrum shows an example where a strong and weak [Nil] and triplets almost 
overlap. The first subtraction is possibly marginally worse than in the previous case. The left spectrum 
shows a case where the lines [Nil] and lines belonging to different triplets are superimposed. The 
automatic procedure makes no attempt to disentangle them. In the first step everything that is at the 
Vmax position of the velocity weight function for [Nil] and Hq, is subtracted. The second subtraction 
finds only a strong [SII] doublet and an Hq, line. 



Fig. 7. — The ip = 90° low resolution spectrum and its reconstruction from velocity components. For 
color coding, abscissa and ordinate range on the reconstruction see Fig. 9. 

Fig. 8. — Trace Sr{X, m) of the spectrum along the line in Fig. la. The [SII] doublet at the red-shift —115 
km/s is unresolved and so is the weak Hq line, but the corresponding [Nil] doublet is wide; another [SII] 
doublet at 814kms~^ (to the right) is clearly about 150kms~^ wide. Note that our spectral resolution of 
16kms~^ is not sufficient to show thermal broadening of or any other line at expected temperatures 
of 7500K. 

Fig. 9. — Reduced low resolution spectra for all slit orientations {ips = 0°, 18° . . . 162°- numbered 
in red). The intensity Wi{p,v), coded according to the color code shown by triangles, is shown as a 
function of velocity (abscissa: v with ticks at lOOkms"^) and position along the slit (ordinate: fh, major 

ticks spaced by 1 arc min). The positive axis of the velocity coordinate points away from the observer 
and the positive axis of the slit coordinate points toward the numbers in Fig. 5. The large orange cross 
indicates the true origin of the coordinate system centered on the pulsar at v = 0, cf. eq. 4. The light 
gray ring with a dot at its center indicates the projection of the spherical shell volume {Rin — 65 arc sec. 
Rout = 90 arc sec) found by the fit. The colors are assigned with respect to the "average bright spot 
mixture of intensities", which is coded gray (center of the triangle, {H : N : S}={0.18 : 0.42 : 0.40}; in 
the average hydrogen is much stronger), while the edges of the triangles correspond to pure H, N or S 
as indicated on top of the figure. Intensity is coded relative to each reduced spectrum as fraction of 
the maximum intensity of the given reduced spectrum. For respective maximum values see text. The 
minimum intensity seen in the electronic version of the image is about 5% of the maximum intensity, 
which is significantly above noise, remember, however, that due to line interference the deconvolution 
is not unique; compare Fig. 7. 

Fig. 10. — Same as Fig. 9 for the high resolution spectra. The S/N for this spectroscopy is quite low (see 
text), so that some specks are due to noise fiuctuations, in particular the trace ip = 18° is the noisiest 
since it has the lowest maximum intensity (see text). This spectroscopy is mainly useful in conjunction 
with low resolution spectroscopy for better velocity calibration. Compare also with Fig. 2. 
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Fig. 11. — A saturated reduced low resolution spectrum for the slit orientation ipg = 144° showing the 
trail of the jet. The color triangle corresponds to the intensity 5%. The signal along the jet is up to 
25% of maximum intensity, yet it is important to realize that the spectrum in this region is crowded, 
so that some hne confusion is probable. Higher resolution and better S/N spectroscopy of this region is 
required. 



Fig. 12. — Crab nebula image reconstruction from the spectral data. Top left: a blend of optical and 
X-ray image, Weisskopf at al. (2000);top ccnter:radio image, Bictcnholz et al. (2001); top right: optical 
image with the dynamical center of the nebula and the projection of the main line emitting lobe; bottom 
left: 6OOOA continuum; bottom right: [Nil] (green), (red), [SII] (blue) spectral components. 



Fig. 13. — An image (a single movie frame) of the Crab nebula from a different perspective. The 
observer on the Earth is in the direction of the —Z axis. The color coding is explained in Sec. 3. 



Fig. 14. — Radial distribution of intensity in the three line sets as a function of the "radial distance" p. 
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